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ABSTRACT 

Kilohcrtz Quasi-Periodic Oscillations(QPOs) in low-mass X-ray binaries (LMXBs) 
may represent the orbital frequencies of innermost Keplerian orbits around accreting 
neutron stars. Attempts have recently been made to derive constraints on the mass 
and the equation of state of the neutron star, by identifying the highest observed 
QPO frequency with the Keplerian frequency at the marginally stable orbit given by 
the general theory of relativity. These estimates have either neglected the effect of 
the neutron star rotation or used an approximate treatment of rotation in general 
relativity. We rederive these constraints using a fully general relativistic formalism 
including the effect of rapid rotation. We also present constraints corresponding to 
the case where the innermost stable orbit touches the stellar surface. 



Key words: Stars: neutron - stars: rotation 
relativity - pulsars: general 



X-rays: stars - equation of state 



9 L,: 
6 



X 



1 INTRODUCTION 

The recent discovery of kilohertz quasi-periodic brightness 
oscillations (kHz QPOs) in certain X-ray binaries (see van 
der Klis 1997 for a recent review) have led to suggestions 
that these are determined by the Keplerian orbital frequency 
near the marginally stable orbit predicted by general rela- 
tivity. Such suggestions have led to attempts to derive con- 
straints on neutron star structure (and possibly also the 
equation of state) using the kHz QPO data (Kaaret, Ford 
& Chen 1997; Zhang, Strohmayer & Swank 1997a; Kluzniak 
1998). It has been assumed in these calculations that the kHz 
QPOs are generated in the innermost stable orbit and that 
the neutron star lies within this orbit. These attempts have 
either neglected the effect of neutron star rotation, or relied 
on the use of an approximate treatment of rotation in gen- 
eral relativity. Since rotation can be important for neutron 
stars in accretion driven old X-ray binaries (Bhattacharya 
& van den Heuvel 1991), a fully general relativistic treat- 
ment of the rotation, in deriving constraints on the neutron 
star structure using kHz QPO data, is imperative. 

In this paper, we explore the consequences of the asso- 
ciation of kHz QPO with Keplerian frequencies at the inner 
edge of the accretion disk, using a fully general relativistic 
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formalism appropriate for a rotational space-time. Unlike 
the previous calculations, we do not make the restrictive as- 
sumption that the radius of the neutron star is less than that 
of the marginally stable orbit. In general, the inner edge of 
the accretion disk need not always be coincident with the 
marginally stable orbit (r m3 ), but can be located anywhere 
outside this radius. If the radius of the neutron star is greater 
than r ms , the innermost possible orbit will be located at the 
surface of the star. Since the structure and r ms for a ro- 
tating neutron star depend on two independent parameters, 
namely, the central density (p c ) and the spin frequency (us), 
a range of values of (p c , us) will exist that will allow solu- 
tions for a Keplerian frequency corresponding to a specific 
kHz frequency. We discuss here two branches of these solu- 
tions: (i) the case where the inner edge of the accretion disk 
is identified with r ms (Solution I) and (ii) the case where 
the disk extends to the surface of the star (Solution II). The 
discussions so far in the literature are confined to cases akin 
to Solution I. This is for two reasons: (a) a constant differ- 
ence between the kHz frequency peaks and its near equality 
to the frequency of the QPO peak observed in X-ray bursts 
are taken to imply that the difference frequency represents 
that of the rotation of the neutron star. The highest QPO 
peak is taken to represent Keplerian frequency at the in- 
nermost marginally stable orbit. Since the frequency of this 
peak is higher than the neutron star rotation frequency, r ms 
is taken to be greater than the radius of the neutron star 
(Kluzniak 1998) and (b) the highest frequency peak is ob- 
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served to be a near constant for different sources and if this 
peak were to originate at the surface of the neutron star, it 
might be difficult to explain the constancy despite the dif- 
ferent magnetic fields and accretion rates that these sources 
may possess (Zhang et al. 1997a). The argument (a) above 
assumes that the accretion disk touches the neutron star 
surface, only if it spins at the centrifugal mass-shed limit 
(in which case, spin frequency w equatorial Keplerian fre- 
quency at the stellar surface). However, since factors such as 
instabilities driven by gravitational radiation reaction (Wag- 
oner 1984; Andersson, Kokkotas & Schutz 1998) may limit 
the spin frequency of the neutron star, there is no compelling 
reason to believe that this assumption is generally valid. The 
arguments in favour of solution I remain source-specific, and 
a main reason for the popularity of this solution is that fea- 
tures relating to the accretion flow and X-ray emission from 
the boundary layer have been computed only for this case 
(Kluzniak and Wilson 1991; Hanawa 1991). Until similar de- 
tailed computations are performed for other situations (e.g. 
solution II), it is not obvious that solution I should always 
be the preferred scenario for kHz QPOs. 



2 CALCULATIONS AND RESULTS 

We look for possible configurations of rotating neutron stars 
in general relativity, whose innermost "allowed" orbits pos- 
sess a QPO frequency of 1220 Hz observed in the source 4U 
1636-53 (Zhang et al. 1997b). This is the maximum value 
of the highest QPO frequency observed so far. Our calcu- 
lations are performed for a representative choice of neutron 
star equation of state (EOS) models, with a view to obtain- 
ing a broad general conclusion. 

There can be two possibilities for the innermost 
marginally stable circular orbits: (i) r ms > R, (ii) r ms < R 
where R is the radius of the neutron star. For the case (ii), 
the innermost "allowed" orbit will be located at the surface 
of the neutron star i.e. r iao = R and for case (i) 
As already mentioned, we designate (i) as Solution I and 
(ii) as Solution II. It may be noted here that the configura- 
tion for which the radius of the innermost stable orbit is just 
equal to the equatorial radius of the neutron star defines the 
boundary between Solution I and Solution II. 

The equilibrium sequences of rapidly rotating neutron 
stars in general relativity and the corresponding values of 
r ms used to derive the results in this paper, have been cal- 
culated using the formalism reported in Datta, Thampan & 
Bombaci (1998). We have chosen the following EOS mod- 
els for neutron star interiors: (i) Bethe & Johnson (1974) 
model V, (ii) Walecka (1974), (iii) Wiringa, Fiks & Fab- 
rocini (1988) model UV14 + UVII (iv) Sahu, Basu & Datta 
(1993) (v) Baldo, Bombaci & Burgio (1997) model BBB2, 
(vi) Bombaci (1995) model BPAL12. In order to conform 
to the notation used by Arnett & Bowers (1977), we refer 
EOS models (i) as (C) and (ii) as (N). The EOS (iii)-(vi) 
are referred to in the text as: (UU), (SBD), (BBB2) and 
(BPAL12) respectively. Of these, (N) and (SBD) are stiff 
EOS, (UU), (C) and (BBB2) are intermediate in stiffness 
and (BPAL12) is a soft EOS. For each of these, we have 
calculated the Keplerian frequency {vk) of a test particle in 
the innermost "allowed" orbit (r iao ) for a range of central 
densities and rotation rates of the neutron star. 
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Figure 1. The Keplerian frequency vk corresponding to the in- 
nermost "allowed" orbit as a function of the gravitational mass 
Mq of the neutron star. The three curves (one solid and two 
dot-dashed) are respectively for three values of neutron star spin 
frequency vg namely 0, 200 and 580 Hz (the last rotation rate 
inferred from 4U 1636-53 by Zhang et al. 1997b). The horizontal 
dashed line corresponds to vk = 1220 Hz, the highest frequency 
QPO observed to date from the X-ray source 4U 1636-53. For all 
portions of the curves above the horizontal dashed line, 1220 Hz is 
obtained at some position exterior to the innermost "allowed" or- 
bit. The values of Mq at intersection points of the curves provide 
the range of masses that admit frequencies of 1220 Hz. 



In Fig. 1 we have illustrated the dependence of the Ke- 
plerian frequency [i>k) on the gravitational mass (Mq) of 
the accreting rotating neutron star for three chosen values 
of its spin frequency (vs): 0, 200 and 580 Hz. In this figure, 
we have considered only three EOS models: (SBD), (UU) 
and (BPAL12), which are respectively, stiff, moderate and 
soft, in order to illustrate the dependence of vk on the EOS. 
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Table 1. Range of masses corresponding to Solution I (intersection of the falling branch of curves in Fig. 1) for various rotation rates 
of the star for vk = 1220 Hz. Only those EOS models which allow Solution I are displayed here. Listed in this table from left to right 
are: the EOS model, central density p c , neutron star spin frequency u$, baryonic mass Mo, gravitational mass Mq, equatorial radius 
R, specific angular momentum j respectively of the neutron star, radius of the innermost marginally stable orbit r ms ht as calculated 
using HT formalism and radius of the innermost "allowed" orbit Ti ao - In this case, Ti ao is equal to the radius of the marginally stable 
orbit dictated by general relativity. Only very slowly rotating configurations (ys < 100 Hz) admit vk = 1220 Hz for EOS model (BBB2) 
and hence there is only one entry against it. The numbers following the letter E in the second column represent powers of ten. 
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The innermost "allowed" stable orbit will be located at the 
neutron star surface for lower values of the mass (the ris- 
ing branch of the curve in Fig. 1 - Solution II) and at the 
marginally stable orbit for higher values of the mass (the 
falling branch in Fig. 1 - Solution I). For a specific value 
of vk (taken here to be equal to 1220 Hz corresponding to 
the highest QPO frequency observed to date in the X-ray 
source 4U 1636-53), Fig. 1 then provides the mass of the 
neutron star (for a given vs and a given EOS): the intersec- 
tion of the curves with the horizontal dashed line defines the 
neutron star configuration which possesses a Keplerian fre- 
quency of 1220 Hz at the innermost "allowed" stable orbit. 
It can be noticed that the peak heights of the curves above 
the horizontal dashed line decreases for increasing stiffness 
of the EOS. It is interesting to see that for model (SBD) 
(a very stiff EOS), the peak value of vk (1259 Hz in this 
case) for v s = 580 Hz is very close to 1220 Hz. Stiff EOS 
models for which this peak lies below vk = 1220 Hz line 
may then be disfavoured. From Fig. 1, further (QPO model 
dependent) constraints on the EOS and neutron star mass 
can be inferred as discussed below: 

(a) If the highest frequency QPOs are produced only at 
the innermost marginally stable orbit (i.e if Solution II is to 
be disfavoured as suggested in Kluzniak 1998; Zhang et al 
1997a; Kaaret et al. 1997), then the constraints obtained on 
the mass can be read off from the higher mass intersection 
points of the curves in Fig. 1. For soft EOS models, such 
an intersection may not obtain. This would then put a con- 
straint against soft EOS models, a point noted in Kluzniak 
(1998). For intermediate and stiff EOS, the inferred mass 
values in this assumption for neutron star spin frequencies 
between and 580 Hz, lies in the range (1.8-2.2) Mq . A 



lower mass limit of 1.8 M using the approximate effects of 
rotation (Hartle-Thorne formalism) was also suggested by 
Kluzniak (1998). If the highest observed frequency in some 
QPO sources have lower values than 1220 Hz, the limits of 
this mass range will increase. For example, for a frequency 
of 1000 Hz, the EOS model (UU) (as well as BPAL12) will 
be disfavoured while the limit obtained from EOS SBD falls 
in the range (2.2-2.5) Mq . 

(b) If the highest kHz QPO peaks are produced only 
at the innermost "allowed" orbit, then another set of con- 
straints: lower limits on the neutron star mass, can be in- 
ferred from Fig. 1. For EOS models (BPAL12) and (UU), 
these limits are in the range (0.9-1.0) Mq and (0.6- 
0.7) M respectively. For EOS SBD, however, this range 
is (1.4-1.8) Mq , which is substantially higher than for the 
other EOS models. For lower values (than 1220 Hz) of the 
highest QPOs, the limits obtained from Solution II would 
allow even lower neutron star masses! 

(c) For the non-restrictive assumption that the highest 
observed QPOs are produced at the inner edge of the ac- 
cretion disk, which might be located outside the innermost 
"allowed" orbit, the only necessary condition imposed will 
be that vk > 1220 Hz at the innermost "allowed" orbit. 
This implies a wider range (from the lower limit of the low 
mass intersection to the higher limit of the higher mass in- 
tersection) of inferred masses. In this case, we cannot rule 
out the validity of the soft EOS models, as the higher mass 
limit will be the maximum mass allowed by the EOS for the 
chosen value of vs- It can also be seen that the constraints 
placed on the mass of the neutron star in this case are not 
as firm as in cases (a) and (b). For EOS models (BPAL12), 
(UU) and (SBD), the mass limits are in the ranges (0.9- 
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Table 2. Range of masses corresponding to Solution II (intersection of the rising branch of curves in Fig.l) for various rotation rates of 
the star for vx = 1220 Hz. The listed quantities are the same as in the first eight columns of Tabic 1. In this case the radius of the star 
is larger than r mB , and Ti ao is located at the stellar surface. 
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1.5) M Q , (0.6-2.1) M and (1.4-2.2) M respectively. For 
lower values of the highest QPO frequency, this range will 
widen making the constraints even less firm, unlike in case 
(a). However, if in any system a QPO frequency is discovered 
which lies above the maximum allowed vk (the peaks in fig. 
1) for the EOS, then the corresponding EOS can certainly 
be ruled out. 

Clearly, the constraints placed on neutron star masses 
(as well as those on the EOS) are dependent on the model 
assumed for the generation of kHz QPOs. 

In Tables 1 and 2 we have listed the values of the neu- 
tron star central density (p c ), spin frequency (fs), the bary- 
onic mass (Mo), the gravitational mass (Mg), the radius 
(R), radius of the innermost allowed orbit (r iao ). For the 



sake of a comparison with the calculations that are based 
on the 'slow' rotation approximation, we have also listed in 
Table 1 r mSi HT and Vk,HT f° r several values of the dimen- 
sionless specific angular momentum j (= J/M G , J being the 
angular momentum of the 'slowly' rotating neutron star). 
The 'slow' rotation approximation is based on the Hartle 
& Thorne (HT) formalism (Hartle & Thorne 1968), which 
treats rotation as a perturbation to the spherically symmet- 
ric space-time. This formalism gives the following analytical 
expressions for the radius of the innermost marginally stable 
orbit and the corresponding Keplerian frequency (Kluzniak 
& Wagoner 1985; Kluzniak, Michelson & Wagoner 1990): 
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fros.HT = 
fK,HT = 
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-3/2 
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(m g ) 
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(2) 



Here c = 1 = G. The HT formalism is valid only for those 
values of the stellar angular velocity that are small in com- 
parison to the centrifugal break-up angular velocity. 

The quantities listed in the tables have been calculated 
by choosing a combination of p c and vs so as to obtain a 
value of vk equal to 1220 Hz. Unlike in Fig. 1, for the results 
that are presented in the tables, no specific value of vs was 
chosen but instead we chose a fixed value for vk = 1220 Hz, 
while exploring various values of vs- Table 1 represents So- 
lution I corresponding to the case Ti ao = r ms and Table 2 
represents Solution II for which ri ao = R- Only those EOS 
models which allow Solution I are displayed in Table 1. For 
EOS (BBB2), Solution I is allowed only for very slowly ro- 
tating configurations with vs < 100 Hz and hence there is 
only one entry against this EOS model in Table 1. For the 
sake of illustration, we include in the Tables 1 and 2 results 
for neutron star rotation frequency of 580 Hz as inferred for 
the source 4U 1636-53 (Zhang et al. 1997b) for all those 
configurations for which such solutions exist. Table 1 shows 
that though r mS: uT varies at most by 10% from ri ao , the dif- 
ferences are substantially larger for the corresponding values 
of vk,ht particularly for large rotation rates. 

It is interesting to ask what the maximum value of vk 
would be for a given EOS model, because this value repre- 
sents the boundary between Solution I and II. In Table 3 we 
list the maximum values of vk along with the corresponding 
Mq and R for different chosen values of vs, namely, 0, 200 
and 580 Hz. It is clear from this table that the softer the 
equation of state, the higher the value of the maximum of 
vk- It may be relevant to note that the EOS model SBD 
has a maximum vk value of 1259 Hz for vs = 580 Hz and a 
corresponding mass of 2.02 Mq . 



3 CONCLUSIONS 

The results of the present study can be summarized as fol- 
lows. We find that the identification of the maximum ob- 
served (so far) QPO frequency (1220 Hz) with vk implies 
the following constraints on the mass of the neutron star 
and the EOS: 

(a) If Solution I is the preferred scenario for kHz QPOs, 
then with the additional assumption that the highest ob- 
served QPO frequencies are produced only at r — r ma , soft 
EOS models for neutron star interiors would be disfavoured. 

(b) If the highest observed QPO frequencies are pro- 
duced at r = ri ao , another set of constraints, possessing 
lower values of mass can be deduced. Together with the 
constraints obtained from case (a) these results limit the 
neutron star mass values to be in either of the two ranges: 
(0.6-1.8) M Q or (1.8-2.2) M . 

(c) If the highest observed QPO frequencies are gen- 
erated at the inner edge of the accretion disk, which may 
be located outside the innermost "allowed" orbit, the con- 
straints implied on the neutron star mass will not be firm 
ones and this value may lie within the range (0.6-2.2) Mq . 
Furthermore, all EOS models allow vk > 1220 Hz values. 



The above results may be considered to be of sufficient 
generality as we have considered EOS models ranging from 
very stiff to very soft. We emphasise that the above con- 
clusions, which constitute a set of constraints, are a con- 
sequence of the inclusion of general relativistic rotational 
effects and the non-restrictive assumption vk > 1220 Hz. 
It is relevant to mention here that the constraints on the 
EOS and the neutron star mass not including the rotational 
effects in general relativity and assuming R < r ms , would be 
somewhat more firm but less general, discounting their util- 
ity. Fig. 1 and Table 3 show that the maximum allowed vk 
for all EOS models that we have considered exceed 1220 Hz. 
We believe, therefore, that determination of the maximum 
QPO frequency and spin rate of the neutron stars among 
a substantially large sample of low-mass X-ray binaries at 
various evolutionary stages in future observations will pro- 
vide more meaningful and general constraints on neutron 
star masses and their EOS. 

Note added: After the submission of this paper, a work on 
a similar theme appeared in the pre-print archives (Miller, 
Lamb & Cook 1998). The results reinforce our conclusions. 
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Table 3. Maximum vk attainable for each EOS for three values of spin frequency vg of the neutron star. Also listed are the corresponding 
values of gravitational mass Mq and radius R. 



EOS u s = u s = 200 u s = 580 

(Hz) (Hz) (Hz) 

uk Mq R ur M~g R uk' Mg R~ 

(Hz) (M Q ) (km) (Hz) (M ) (km) (Hz) (Mq ) (km) 



(C) 


1783 


1.224 


10.898 


1819 


1.284 


10.884 


1817 


1.372 


11.122 


(N) 


1560 


1.330 


12.250 


1606 


1.468 


12.364 


1584 


1.619 


12.857 


(UU) 


1743 


1.244 


11.122 


1774 


1.320 


11.170 


1769 


1.439 


11.497 


(SBD) 


1303 


1.692 


15.008 


1306 


1.800 


15.208 


1259 


2.021 


16.235 


(BBB2) 


1732 


1.258 


11.218 


1760 


1.317 


11.226 


1770 


1.436 


11.494 


(BPAL12) 


1752 


1.254 


11.106 


1800 


1.293 


10.989 


1840 


1.373 


10.974 
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